Frenkel-Kontorova model with cold trapped ions 
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We study analytically and numerically the properties of one- dimensional chain of cold ions placed 
in a periodic potential of optical lattice and global harmonic potential of a trap. In close similarity 
with the Frenkel-Kontorova model, a transition from sliding to pinned phase takes place with the 
increase of the optical lattice potential for the density of ions incommensurate with the lattice 
period. Quantum fluctuations lead to a quantum phase transition and melting of pinned instanton 
glass phase at large values of dimensional Planck constant. The obtained results are also relevant 
for a Wigner crystal placed in a periodic potential. 
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Nowadays experimental techniques allow to store thou- 
sands of cold ions and observe various ordered structures 
formed by Coulomb repulsion in ion traps yj. These 
structures include the one-dimensional (ID) Wigner crys- 
tal, zig-zag and helical structures in three dimensions. 
The Cirac-Zollcr proposal of quantum computations with 
cold trapped ions generated an enormous experimen- 
tal progress in this field with implementations of quan- 
tum algorithms and quantum state preparation with up 
to 8 qubits 0. [l| . In these experiments 0, ions form 
a ID chain placed in a global harmonic potential which 
frequency u> determines the eigenfrequencies of chain os- 
cillations being independent of ion charge 0, |(| . Highly 
accurate experimental control of the chain modes allows 
to perform quantum gates between internal levels of ions. 
In addition to ion traps, modern laser techniques allow 
to create periodic optical lattices and store in them thou- 
sands of cold atoms (see e.g. 0). A single ion dynamics 
in an optical lattice has been also studied experimentally 
8] . The combination of these two techniques makes pos- 
sible to study experimentally the properties of a ID chain 
of few tens of ions placed in an optical lattice and a global 
harmonic potential at ultra low temperatures. 

In this Letter we analyze the physical properties of 
such a system and show that it is closely related to 
the Frenkel-Kontorova Model (FKM) [£j which gives 
a mathematical description of various physical phe- 
nomena including crystal dislocations, commensurate- 
incommensurate phase transitions, epitaxial monolayers 
on a crystal surface, magnetic chains and fluxon dynam- 
ics in Josephson junctions (see [lOj and Refs. therein). 
As in the classical FKM the ion chain exhibits the Aubry 
analyticity breaking transition when the amplitude 
of optical potential becomes larger than a certain crit- 
ical value. Below the critical point the classical chain 
can slide (oscillate) in the incommensurate optical lat- 
tice while above the transition it becomes pinned by the 
lattice and a large gap opens in the spectrum of phonon 
excitations. Above the transition the positions of ions 



form a devil's staircase corresponding to a fractal Can- 
tor set which replaces a continuous Kolmogorov-Arnold- 
Moser (KAM) curve in the phase space below the transi- 
tion. According to yjj the FKM ground state is unique 
but in the pinned phase there are exponentially many 
stable configurations which are exponentially close to the 
ground state and the system resembles a fractal spin glass 
In the quantum FKM [HQ El 

the instanton tun- 
neling between quasidegenerate configurations leads to a 
melting of the pinned instanton glass and the quantum 
phase transition into sliding phonon phase takes place 
above a critical value of the effective Planck constant 
h 0| . In contrast to the FKM with nearest neighbor 
interactions between particles, the ions have long range 
Coulomb interactions. In spite of that we show that the 
ion system is effectively described by the FKM and has 
similar phase transitions which can be studied experi- 
mentally. To analyze the properties of the FK Ion Model 
(FKIM) with few hundreds of ions we use various nu- 
merical methods including classical and Quantum Monte 
Carlo (QMC) algorithms as described in [H [3 El- 
The dimcnsionless FKIM Hamiltonian has the form: 
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where Pi,Xi are ion momentum and position, K gives 
the strength of optical lattice potential and all N ions 
are placed in a harmonic potential with frequency to. To 
make a transfer from to dimensional physical units 
one should note that the lattice constant d in K cos(xi/d) 
is taken to be unity, the energy E — H is measured 



in units of ion charge energy e jd and uj 



2 d 3 /e 



where m is ion mass. In the quantum case = —ihd/dxi 
with dimensionless h measured in units H — > H/(ey/md). 
In the quantum case we use the approximation of distin- 
guishable ions which is well justified when the distance 
between ions imposed by the harmonic potential is com- 
parable with the lattice period @|. 

We start the discussion from the classical case. Here 
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FIG. 1: (Color online) Functions related to the dynamical 
map © obtained from the ground state equilibrium positions 
xi of N = 50 ions for u = 0.014 (lo 2 = 1.986 x 10~ 4 ) at K = 
0.03 (open circles, left column) and K — 0.2 (full circles, right 
column). Panels show: the kick g(x) function (top); the phase 
space (p,x) of the map @ with g(x) = — sin a; (green/gray 
points) and actual ion positions (red/black circles) (middle); 
the hull function h(x) (bottom). The ion positions are shown 
as x — Xi(mod2n) for the central 1/3 part of the chain. 



the stable configurations with minimal energy have Pj = 
and satisfy the conditions dH/dx.i = 0. In approx- 
imation of only nearest neighbor interacting ions these 
conditions lead to the dynamical recursive map for equi- 
librium ion positions xc 



Pi+i = Pi + Kg(xi) , Xi 



where the effective momentum conjugated to Xi is pi — 
l/(x{ — Xi_i) 2 and the kick function Kg(x) = —lo 2 x — 
K sin a;. To check the validity of this description we find 
the ground state configuration using numerical methods 
described in 0, ^| . The harmonic frequency lo is cho- 
sen in such a way that in the middle of the chain the ion 
density v = 2ir/ (x\ — xq) at K — is equal to the golden 
mean value v = v g = (V5 + l)/2. This corresponds to 
an incommensurate phase with the golden KAM curve 
usually studied for the Aubry transition 0, [n| . For a 
fixed v ~ 1 the strength of harmonic potential lo 2 drops 
with N according to relations found in 0, Q . The nu- 
merical values Xi allow to determine the function g{x) 
using relations J2J and show that with good accuracy 
g(x) = — s'mx (see Fig^). We attribute this result to an 
effective screening of global harmonic potential by long 
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FIG. 2: (Color online) Phonon spectrum ui(k) as a function of 
scaled mode number k = i/N (i = 0, . . . , N — 1), for K — 0.03 
(bottom curve, squares) and K — 0.2 (top curve, points) for 
the case of FigQ 
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FIG. 3: (Color online) Minimal excitation frequency u>o(K) 
as a function of periodic potential strength K for the golden 
mean ion density v g and number of ions N — 50 (squares; 
U) = 0.014, = 1.986 x 10" 4 ), N = 150 (circles; u = 
0.00528, a; 2 = 2.794 x 10" 5 ), N = 300 (triangles, uo = 
0.00281, a; 2 = 7.913 x 10~ 6 ). The critical point K c ~ 0.05 



(2) is marked by arrow; inset shows data near K c 



range interactions between ions in the central part of 
the chain where the ion density is approximately con- 
stant. For small values of optical potential (K = 0.03) 
the chain is in the sliding phase with a continuous KAM 
curve in the plane (p, x) while above the Aubry transi- 
tion at K = 0.2 it is in the pinned phase and forms a 
fractal Cantor set in the phase space (FigQ. The qual- 
itative change of chain properties is also seen via the 
hull function h(x) which gives the ion positions in peri- 
odic potential vs. unperturbed positions at K = both 
taken mod 2tt: for K < K c one has a continuous func- 
tion h(x) rs x while for K > K c the hull function has a 
form of devil's staircase with clustering of ion positions 
at certain values (Fignj. More data for a larger number 
of ions N = 150 are shown in Appendix (see Fig@. 
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FIG. 4: (Color online) Number of equilibrium configurations 
N s as a function of their relative excitation energy AU/U 
above the ground state for 50 (blue/black) and 150 (ma- 
genta/gray) ions at K = 0.2 with ui as in Fig[3] (logarithms 
are decimal). 




The structural changes in the chain below and above 
the Aubry transition at the critical value K c are also 
clearly seen in the spectrum of phonon excitations shown 
in Fig|3 (see also Fig0in Appendix for N = 150). For 
K < K c the spectrum of phonon modes has a sound 
like form at small wave vectors fc, going down to a min- 
imal oscillation frequency loq ss lo, which is small com- 
pared to frequencies of periodic potential. For K > K c 
the spectrum of excitations is characterized by a phonon 
gap when loo becomes independent of number of ions N 
(FigOJ. According to the data of FigOU for loq(K) and 
of FiglHI in Appendix for the hull function the critical 
point is located at K c w 0.05. The value K c can be 
obtained approximately by a reduction of map @ to 
the Chirikov standard map (3- F° r that the second 
equation can be linearized near l/yf> = 27r/V ff that 
gives the Chirikov standard map with chaos parameter 
K eff = K(2ir/v g Y/2 and the critical value K c = 0.034 
corresponding to K e ff — 1. This value is smaller than 
the value K c rs 0.05 obtained from numerical data that 
can be attributed to a strong dependence of K e ff on v g . 

The obtained results show that the ion chain in a pe- 
riodic potential has close similarities with the FKM. The 
map J2J) describes quite accurately the central 1/3 part 
of the chain where the ion density is approximately con- 
stant. Outside of this part the ion density starts to grow 
and deviations from the map start to be visible in (p, x) 
plain and h(x), even if the ^-function remains rather close 
to g = — sin x. A separation on two parts is also seen 
in the lowest phonon eigenmodes which are localized in 
the central part of the chain (see Figs 191101 of Appendix). 
Also, in similarity with the FKM [12| . we find that in the 
pinned phase at K > K c the FKIM has properties of a 
spin glass with enormous number of stable equilibrium 
configurations N s in the very close vicinity of the ground 



FIG. 5: (Color online) Formfactor F(k) (see text) of the 



chain with TV = 50 ions and lo — 0.014, lo = 1. 
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(a): The classical incommensurate phase at K — 0.03, fi — 0, 
arrows mark the peaks at integer multiples of golden mean 
density v g . (b) The pinned phase case at K = 0.2 for h — 
(bottom black curve), ft = 0.1 (middle red curve shifted 10 
units upward), h — 2 (top green curve shifted 20 units upward, 
for clarity F(k) is multiplied by factor 5). The temperature of 
the quantum chain is T — h/ro with to = 400 so that T <C K 
and T < hu) {K). 



state (see Fig0J. 

These quasidegenerate configurations become very im- 
portant in the quantum chain where quantum tunnel- 
ing between them leads to nontrivial instanton excita- 
tions. To study the quantum case we use the QMC 
approach Iltl Hi} with the Euclidean time r in 
the interval [0, To] and system temperature T = K/tq. 
We use large values of r (e.g. t > 400) so that 
the temperature T is small compared to periodic po- 
tential K and phonon gap Hloq. To see the structural 
changes produced by quantum fluctuations it is conve- 
nient to compute the formfactor of ion positions defined 

l 1 2 

asF(fc) = | exp(ikXi(T)) | )/SN, where summation is 
taken over the central part of the chain with 6N w NJ 3 
ions, Xi(r) are ions positions at Euclidian time r and 
brackets mark averaging over r. In the classical case 
Xi are the equilibrium positions of ions. For h = the 
formfactor peaks are equidistant in v g that clearly show 
the incommensurate sliding phase at K = 0.03 < K c 
(FigEt). For K = 0.2 > K c the peaks in F(k) are at 
integer k values showing that ions are pinned by the op- 
tical lattice (FiglSb bottom). In the quantum case at 
K = 0.2 > K c and small values of h = 0.1 the density of 
instantons is small and quantum tunneling does not de- 
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stroy the pinned phase. However, at large h = 2 quantum 
fluctuations lead to melting of pinned phase, the integer 
peaks disappear and F(k) becomes continuous (FigEb). 
This means that the ion chain can slide in the optical 
lattice. The quantum phase transition from pinned in- 
stanton glass to sliding phonon gas takes place at zero 
temperature and certain K c ~ 1. 

The formfactor F(k) can be experimentally obtained 
by light scattering on the ion chain and in this way vari- 
ous phases of the system can be detected. For experimen- 
tal conditions like in 0] the distance between Ca ions is 
about Ax « 5/itm so that the golden mean density v g cor- 
responds to the lattice constant d — v g Ax/2Tr w 1.3/zm. 
This value can be realized by laser beams crossed at a 
fixed angle. The transition to pinned phase takes place at 
the optical lattice potential V = K c e 2 /d « 0.6K that in 
principle can be reached in strong laser fields. For 40 Ca + 
ions with such d the dimensionless Planck constant is 
% e jf = hj (ey/md) ~ 3 x 1CP 5 . This means that such ex- 



periments with cold trapped ions can be performed in a 
deep semiclassical regime hardly accessible to the QMG 
numerical simulations. We also note that a strong gap in 
phonon spectrum u>k for K > K c (Figs. EH may be use- 
ful for protection of quantum gate operations in quantum 
computations. 

Higher values of dimensionless h can be reached with 
electrons forming a Wigner crystal placed in a periodic 
potential where h « 0.1 for d m bum. Such a situation 
may appear in ID electron wires, molecular structures 
or nanotubes where a devil's staircase behavior has been 
discussed recently [lij . In the regime when the number 
of electrons per period is of the order of v ~ 1 the effects 
related to statistics of particles are not so crucial and the 
quantum melting of the pinned phase should qualitatively 
follow the scenario described here. Hence, the obtained 
results describe also a more general problem of a Wigner 
crystal in a periodic potential. 

This work was supported in part by the EC IST-FET 
project EuroSQIP. 
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FIG. 6: (Color online) Same as in FigQbut for N — 150 ions 
at lu = 0.00528 {lu 2 = 2.794 x 10" 5 ) corresponding to the 
golden mean density of ions in the middle of the chain. 




APPENDIX 



Here in Figs 161 ill we present additional data showing 
results with larger number of ions. 



FIG. 7: (Color online) Phonon spectrum u)(k) as a function of 
scaled mode number k = i/N (i = 0, . . . , N — 1), for K — 0.03 
(bottom curve, squares) and K = 0.2 (top curve, points) for 
the case of Fig|S]at N = 150 (compare with data of FigElfor 
N = 50). 



* |http : //www.quantware .ups-t lse . f r| 
[1] G. Birkl, S. Kassner and H. Walther, Nature 357, 310 
(1992). 

[2] J.I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 

[3] H. Haffner, W. Hansel, C.F. Roos, J. Benhelm, D. Chek- 
al-kar, M. Chwalla, T.Korber, U.D. Rapol, M. Riebe, 
P.O. Schmidt, C. Becher, O. Ghne, W. Diir and R. Blatt, 
Nature 438, 643 (2005). 



[4] D. Leibfried, E. Knill, S. Seidelin, J. Britton, 
R.B. Blakestad, J. Chiaverini, D.B. Hume, W.M. Itano, 
J.D. Jost, C. Langer, R. Ozeri, R. Reichle and 
D.J. Wineland, Nature 438, 639 (2005). 

[5] D.H.E. Dubin and T.M. O'Neil, Rev. Mod. Phys. 71, 87 

(1999) . 

[6] G. Morigi and S. Fishman, Phys. Rev. Lett. 93, 170602 
(2004). 

[7] I. Bloch, T.W. Hansen and T. Esslinger, Nature 403, 166 

(2000) . 



5 



it 
i/2 
X o 

-ji/2 



it 



Jt/2 

-Jt/2 
-it 



(b) 



*** 



(d) 



-it/2 7t/2 

X; mod(2ji) 



it -it -it/2 it/2 it 

X; mod(2ji;) 



FIG. 8: (Color online) Evidence of the Aubry transition in 
the hull function. The hull function h(x) shows the ion po- 
sitions Xi(K) mod (27r) in presence of periodic potential as 
a function of ion positions x[ ^ mod(27r)) without potential 
(K = 0). Here N = 150, uj = 0.000528, uj 2 = 2.794 10~ 5 , (a) 
K = 0.04; (b) K = 0.055; (c) K = 0.065; (d) K = 0.1. Only 
the central 1/3 part of the chain is shown. 
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FIG. 9: Phonon eigenmodes in the chain with N — 50 
ions at the golden mean ion density with uj = 0.014, a; 2 = 
1.986 x 10~ 4 . Panels show the amplitude of eigenmode (in ar- 
bitrary linear units) vs. ion positions Xi varied from minimal 
to maximal value for lowest 4 modes with number 0, 1, 2, 3 
(from left to right). Top row is for K = 0.03 and bottom row 
is for K = 0.2. 
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FIG. 10: Same as in FigHfor = 150, to = 0.00528, 
2.794 x 10" 5 . 
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FIG. 11: (Color online) Formfactor F(k) (see text) of the 
chain with TV = 150 ions and u = 0.00528, lo 2 = 2.794 x 10" 5 . 
(a): The classical incommensurate phase at K — 0.03, h — 0, 
arrows mark the peaks at integer multiples of golden mean 
density v g . (b) The pinned phase at K — 0.2 for h = 
(bottom black curve), h = 0.1 (middle red curve shifted 20 
units upward), fi — 2 (top green curve shifted 40 units upward, 
for clarity F(k) is multiplied by factor 5). The temperature of 
the quantum chain is T — ft/ro with to = 400 so that T <C K 
and T < hu) (K). 
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